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ABSTRACT 

The Thermal Control Subsystem (TCS) is a key part of the Advanced Plant Habitat (APH) for the 
International Space Station (ISS). The purpose of this subsystem is to provide thermal control, 
mainly cooling, to the other APH subsystems. One of these subsystems, the Environmental 
Control Subsystem (ECS), controls the temperature and humidity of the growth chamber (GC) 
air to optimize the growth of plants in the habitat. The TCS provides thermal control to the ECS 
with three cold plates, which use Thermoelectric Coolers (TECs) to heat or cool water as needed 
to control the air temperature in the ECS system. In order to optimize the TCS design, pressure 
drop and heat transfer analyses were needed. 

The analysis for this system was performed in Siemens NX Thermal/Flow software (Version 8.5). 
NX Thermal/Flow has the ability to perform ID or 3D flow solutions. The ID flow solver can be 
used to represent simple geometries, such as pipes and tubes. The ID flow method also has the 
ability to simulate either fluid only or fluid and wall regions. The 3D flow solver is similar to 
other Computational Fluid Dynamic (CFD) software. 

TCS performance was analyzed using both the ID and 3D solvers. Each method produced 
different results, which will be evaluated and discussed. 


INTRODUCTION 

Advanced Plant Habitat is a large, closed loop plant growth chamber that is designed to be used 
in an EXPRESS rack for the International Space Station (ISS). Environment within the growth 
chamber (GC) will be controlled, such as temperature and humidity, using the Environmental 
Control System (ECS). Finally thermal control is provided to the ECS using the Thermal Control 
Subsystem (TCS), which is a series of cold plates and Thermoelectric Coolers (TECs) working in 
unison with a pressurized water loop. 

The TCS is broken into three different components. The first is the thermal control unit initial 
(TCUI). This is made from 6061-T6 aluminum and has eight thermoelectrics attached to one side 
of the device. Each TECs is set to generate 18 watts of heat, used to cool the air and heat the 
fluid. 

The second component is the humidity control unit (HCU). This is made from 6061-T6 aluminum 
and has 16 TECs. The TECs are attached to one side of the device, each generating 4.9 watts of 
heat, used to cool the air and heat the fluid. The final component of the TCS is the thermal 
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control unit final (TCUF). This is also made from 6061-T6 aluminum and has eight TECs attached 
to the top of the unit. Each TEC on this unit is set to generate 6.1 watts of heat, used to heat the 
air and cool the fluid (opposite of the TCUI and HCU). A combined image of the TCS system, 
which includes all three components is displayed in Figure 1. 



Figure 1. Image of Thermal Control Subsystem 


DETAILED ANALYSIS 

The analysis completed on the TCS was done using two different methods in Siemens NX 
Advanced Thermal/Flow (Version 8.5). The first method used only NX Advanced Thermal. This 
modeled the cold plates as 3D tetrahedral elements, while simplifying the fluid region and tube 
wall to a ID mesh (ID flow solver). The second methods used NX Advanced Thermal & Flow and 
modeled the cold plates, tube wall, and fluid region using 3D tetrahedral elements (3D flow 
solver). 


Method 1: ID flow solver 

Using this method the cold plates were meshed in 3D tetrahedral elements, while the fluid 
region and tube wall were meshed using a ID mesh. The thermal and fluid regions were 
connected using convection coupling methods. This couple's convection between the polygon 
surfaces of the cold-plates, to the curves of the ID fluid region, also called the duct flow model. 
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Additionally the boundary conditions and solver inputs to this method included the following: 


• Constant mass flow rate for each flow case 

• Serial processor solving 

• Constant heat load for each TEC 

• Rough walled tubing 

• Turbulence begins at Reynolds > 2400 (NX default) 

• Forced convection coupling 

• Steady state solutions 



Figure 2. TCS Mesh using Method 1 (ID flow solver) 

Method 2: 3D flow solver 

In this method 3D tetrahedral elements were used to mesh the cold plates, tube wall, and fluid 
regions. For the fluid region a structured quadrilateral o-grid would typically be used to resolve 
the boundary layers. Instead an unstructured tetrahedral mesh was used in conjunction with 
wall functions. Using wall functions allows for the mesh to be larger in size, so the viscous 
sublayer does not need to be resolved. If wall functions were not used, a quadrilateral mesh 
would have been implemented and a y+= 1 would have been desired. An image of the tube 
wall and fluid region mesh is shown in Figure 3. 
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Figure 3. TCS Fluid Region Mesh using Method 2 (3D flow solver) 

The y + was calculated using [1] estimated for each turbulent flow case, which was assumed to 
be when the Reynolds number was greater than 2000. The equation typically used to 
determine this may be seen in Equation 1. 


y + = 


y-u* 

v 


Equation 1 


From [2] pages 358-359 

The purpose of calculating the y + values for these flow cases was to determine the mesh size 
that would be required. An image of the turbulent boundary layer from [3] is shown below in 
Figure 4. In this figure we are concerned with the inner layer shown by the #4. This is compose 
of the viscous sublayer #1 (y + < 5), the buffer layer #2 (5 < y + < 30), and the log-law region 
#3 (30 < y + < 200). When using wall functions you must be sure that the y + produced by 
the model is greater than 5 and in the buffer layer or higher. The main goal of using wall 
functions is to be outside the viscous sublayer. For this analysis hand calculations were initially 
performed to try and keep the mesh size greater than a y + of 30 using the average velocity. 
After calculations were performed, it was determined that the y + value of the model ranged 
between 7.8 to 13.5, which was lower than desired, but still acceptable. The y + values 
determined at the model inlet are shown in Figure 5. 

Additionally the boundary conditions and solver inputs to the model included the following: 

• Constant mass flow rate for each flow case 
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• Parallel processor solving (6 processors) 

• Mixing Length Turbulent model (Reynolds >2000) 

• Constant heat load for each TEC 

• Rough walled tubing 

• Mesh mating conditions 

• Steady state solutions 



Figure 4. Turbulent Boundary Layer Depiction from [3] 



Figure 5. Y+ Values for inlet of TCUI for Method 2 
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RESULTS 


Originally the two different methods were intended to be used for results verification. After 
initial results were completed, the results for cold plate temperature and fluid pressure drop 
were scene to have large discrepancies, while the fluid temperature exiting the TCS was not 
seen to change. The following sections will give an overview of the results determined from the 
analysis. 


Method 1: ID flow solver results 

The exact results were extracted from contour profiles shown in the Siemens NX software. The 
average temperature on all three cold plates may be seen in Figure 6. This shows that the 
average temperatures of the TCUI and HCU start high and begin to reduce as flow rate 
increases, while the TCUF has the opposite response (due to TECs cooling). While the trending 
of the results looks reasonable, the magnitude of the cold plate temperature difference 
seemed large and will be compared to Method 2. The temperature contour profiles at 70 
Ibm/hr may be seen for all three components of the TCS in Figure 7. 

In addition to the thermal surface temperatures of the cold plates, the fluid exit temperature 
and fluid pressure drop were also evaluated. The fluid temperature exiting each cold plate may 
be seen in Figure 8, while the fluid pressure drop through the TCS is shown in Figure 9. Both 
figures appear to have realistic trending results, however when the pressure drop numbers are 
compared to Method 2 a larger discrepancy exists. The total pressure drop across the TCS and 
several other subsystems must remain below 2.8 +/- 0.15 PSID, which is the system 
requirement. Pressure drop results from this method appear to be well within that range, 
showing a maximum pressure differential of 0.243 PSID. These results will be compared against 
those received from 3D flow model (CFD) used in Method 2. 
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tcu full 021814 asm slmlC : MF 70 Ibm-hr Result 
Load Case 1 , Static Step 1 
Temperature - Elemental, Scalar 
Min : 65.49, Max : Ujy^LUnits * F 


tcu full 021814 asm simIC : MF 70 bnvhr Result 
Load Case 1. State Step 1 
Temperature - Elemental, Scalar 
Min : 65.49, Max 1 16.33, Units = F 



Figure 7. Temperature Contour of TCS Cold Plates for Method 1 @ 70 Ibm/hr 
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Figure 8. Fluid Temperature Exiting TCS System for Method 1 
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Figure 9. Fluid Pressure Drop in TCS for Method 1 
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Method 2: 3D flow solver results 


The exact results were extracted from contour profiles shown in the Siemens NX software for 
both 3D thermal and fluid meshes. The average plate surface temperatures for all three cold 
plates have been plotted in Figure 10. When comparing these results with those shown 
previously in Figure 6 some discrepancies are noticed. Both Figure 6 (Method 1) and Figure 10 
(Method 2) have surfaces temperatures decreasing as fluid travels from the TCUI to TCUF. Both 
methods use the same TECs heat input and cold plate mesh. The thermal and fluid models for 
Method 2 are connected with mesh mating conditions, instead of convection coupling. This 
could attribute to the different results. The temperature contour profiles at 70 Ibm/hr may be 
seen for all three components of the TCS in Figure 11. 

As with Method 1, the fluid exit temperature and fluid pressure drop through each cold plate 
was captured in Figure 12 and Figure 13 for Method 2. By comparing the fluid exit temperature 
in Figure 8 (Method 1) and Figure 12 (Method 2) shows very little difference between the two 
graphs. Results appear to be similar trending and have similar temperature values. Next the 
pressure drop differences were compared between Figure 9 (Method 1) and Figure 13 (Method 
2). Results from Method 2 have a larger overall pressure drop then in Method 1. The increased 
pressure drop from Method 2 to 1 is about a 246% increase. Furthermore cases for Method 2 
(specifically at 80, 90, and 100 Ibm/hr) were turbulent, and turbulence modeling was used. The 
two models tested were the Mixing Length (single equation model) and the K-e (two-equation 
model). The variation in pressure drop using these models was found to be less than 3% and 
the higher pressure drop results were displayed here (Mixing Length model). 

While the fluid temperatures in both analyses are almost identical, other parameters differ 
greatly. Pressure drop in the fluid is higher with the 3D modeling, which is expected with 
turbulent energy is being modeled. Furthermore there is a very large difference in the 
predicted plate temperatures. These different methods show two very different temperature 
profile trends. This difference could be related to issues coupling fluid and thermal interactions 
and these models should be validated with acceptance test data. 
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Figure 11. Temperature Contour of TCS Cold Plates for Method 2 @ 70 Ibm/hr 
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Fluid Temperature 



— TCUI-Out 
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Figure 12. Fluid Temperature Exiting TCS System for Method 2 



Figure 13. Fluid Pressure Drop in TCS for Method 2 
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TEST DATA COMPARISON 

Currently integrated system testing is being completed at Orbital Technologies (Orbitec) in 
Wisconsin. While no final results have been published, some preliminary testing results have 
been provided to the author for model comparison purposes. Based on data provided the 3D 
CFD model used in Method was closer to experimental data. The difference in pressure drop 
results between the CFD model and experimental data ranged between 4% to 50% depending 
on flow rates, while percent difference between the ID duct network method and 
experimental data were much higher. 


CONCLUSIONS 

The results from the two different modeling methods (ID vs. 3D flow model) produced similar 
fluid exit temperature results, while the cold plate temperatures and fluid pressure drop results 
were very different. After reviewing system requirements it was determined that a higher 
importance should be placed on the pressure drop results, instead of temperature. The TCS is 
only a single subsystem in a larger system that all must remain below a total pressure drop of 
2.8 +/- 0.15 PSID. Results from the 3D flow model showed the TCS using up to 21.4% of the 
allowable pressure loss, while the ID flow model used only 8.4%. For these reasons the 3D flow 
solver results were used for developing the prototype TCS for testing. 

Since the completion of this analysis, integrated testing has been performed on the TCS at 
Orbital Technologies (Orbitec) in Wisconsin. Preliminary testing results showed that the results 
from the 3D flow model were similar to actual system performance. Release of the results from 
Orbitec are still pending. 


NOMENCLATURE 

APH advanced plant habitat 

EXPRESS Expedite the Processing of Experiments for Space Station Racks 
°F degree Fahrenheit 

ISS International space station 

Ibm/hr pounds of mass per hour 

PCG plant growth chamber 

PSID pound per square inch differential 

TCS thermal control system 

TEC thermoelectric cooler 

TCUF thermal control unit final 

TCUI thermal control unit initial 

HCU humidity control unit 

y distance to the nearest wall 

y + Non-dimensional wall distance 

u* friction velocity 

V kinematic viscosity 
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Background 

• Thermal Control System (TCS) provides thermal control, 
mainly cooling, to several Advanced Plant Habitat (APH) 
subsystems. One of these systems is the Environmental 
Control System (ECS), which provides temperature and 
humidity control inside the growth chamber to optimize 
plant growth. 

• TCS consists of three cold plates that use 
Thermoelectric Coolers (TECs) to heat or cool a water 
loop as needed to meet ECS temperature requirements. 
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Thermal Control System (TCS) 




Humidity Control Unit 
(HCU) 


Thermal Control Unit Final 
(TCUF) 
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Problem Statement 


• Pressure drop and heat transfer analysis was performed 
to optimize the TCS design. 

• Computational Fluid Dynamics (CFD) was performed in 
Siemens NX Thermal/Flow. 

- NX Thermal (Only): Includes 3D mesh thermal solver and ID 
mesh solver for fluid flow in pipe or tubes. This may include wall 
thickness (if desired). 

- NX Thermal/Flow: Includes 3D mesh thermal solver and 3D CFD 
flow solver. 

- Both methods were used and the results will be compared in this 
presentation. 
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Model Differences 




1 D flow model 

• Steady state solution 

• Constant mass flow rate 

• Constant heat loads (per 
TEC) 

• Rough walled tubing 

• Turbulence begins at 
Reynolds > 2400 (default) 

• Forced convection 
coupling 

• Single processor 


3D flow model 

• Steady state solution 

• Constant mass flow rate 

• Constant heat loads (per 
TEC) 

• Rough walled tubing 

• Turbulence begins 
Reynolds > 2000 

• Mesh mating conditions 

• Uses wall functions 

• Parallel (6 processors) 
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Method 1: ID duct flow network analysis 







Models fluid region as one- 
dimensional fluid duct. 

- Pressure 

- Velocity 

- Mass flow rate 

- Temperature 

1 D mesh 

- Fluid only 

- Fluid and tube wall 

Heat transfer (Thermal coupling) 

- Conduction, convection, radiation 

- Free, forced and user defined 
convection 
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Method 1: ID duct flow network analysis 


Does not use NX Flow 

- Only NX Thermal is used 

Duct boundary conditions 

- Fan/pump 

- Opening 

- Pressure 

- Flow properties 

- Duct to 3D flow connection 

• Connects 1 D duct to 3D CFD model 



TFAWS 2014 -August 4-8, 2014 




Method 2: 3D CFD Analysis 





3D mesh 

- Structured hexahedral 

• Quad 2D face and sweep (manual) 

- Unstructured tetrahedral 

• Either automatic or manual meshing 

Turbulent model 

- K-Omega 

- K-Epsilon 

- Mixing Length 

- LES 

Both NX Flow and Thermal 
Serial & Parallel processing 
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Method 2: 3D CFD Analysis 




Parallel solver used fully 
coupled pressure velocity 
scheme 

All results converged for a 
steady state solution 

Residuals used for the analysis 
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Method 1: Duct Meshing Sizing 

• 1 D duct network included both tube wall and fluid 

- Fluid and 3D thermal model for wall were tested 

• Tube used 0.075 inches element size 

• ID duct mesh (fluid) 

- 0.125 inches element size 

• 1 D duct mesh (tube wall) 

- 0.125 inches element size 

• 3D tetrahedral cold plate mesh 

- 0.125 inches element size (0.395 inch recommended) 

• Mesh sensitivity was tested at half and double of the 
element size and found the changing mesh size had no 
effect on the solution 
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Method 1 : Duct Mesh 




'i 
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Method 1 : Duct Mesh 



ID duct with wall 


ID duct with 3D wall 
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Method 2: CFD Mesh Sizing 


Fluid mesh (3D tetrahedral) 

- 0.33 inch element size 

- wall functions were used so viscous sublayer did not need to be 
resolved 

Tube wall mesh (3D tetrahedral) 

- 0.075 inch element size (0.164 inch recommended) 

Cold plate mesh (3D tetrahedral) 

- 0.125 inches element size (0.395 inch recommended) 

Mesh sensitivity was tested by varying the mesh size 
and no effects to the solution was noticed 
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Inner layer 

1 . Viscous sublayer 

2. Buffer layer 

3. Log-law region 

Goal was to keep 
mesh large enough 
so y+>30 

- Actual ranged 7.8 to 
13.5 


TFAWS 2014 


CFD Mesh 




August 4-8, 2014 


14 



3D flow model mesh 
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Cold Plate Surface Temperature Results 

1 D flow model 3D flow model 
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Fluid Temperature Results 



1 D flow model 


3D flow model 
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Fluid Pressure Drop Results 



1 D flow model 


3D flow model 
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Conclusions 


• Models produced similar fluid temperatures, but differed 
in pressure drop and surface temperature results. 

• Testing sensitivity in Method 1 to convection coupling or 
Method 2 to turbulence model did not produce drastic 
changes 

- Grid size sensitivity also did not effect the solution 

- Additional pressure drop could be do to additional swirling 
velocity terms (problems would be in any ID network model) 

• Preliminary data from the integrated testing at Orbitec 
showed the CFD method was closer to in actual 
pressure drop to prototype. 

- 1 for 1 comparisons aren’t available since on the spot changes 
were made to the prototype were made, moving the 
configuration away from the geometry analyzed. 
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Questions? 
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Backup 



1 D flow model (wall included) 


1 D flow model (wall excluded) 
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Backup Flow Model Meshes 




1 D flow model 3D flow model 
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Backup 



1 D flow model 


3D flow model 
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